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We consider zero temperature behavior of dynamic response functions of ID systems near edges of 
support in momentum-energy plane {k,uj). The description of the singularities of dynamic response 
functions near an edge e{k) is given by the effective Hamiltonian of a mobile impurity moving in 
a Luttinger liquid. For Galilean-invariant systems, we relate the parameters of such an effective 
Hamiltonian to the properties of the function e{k). This allows us to express the exponents which 
characterize singular response functions of spinless bosonic or fermionic liquids in terms of e{k) and 
Luttinger liquid parameters for any k. For an antiferromagnetic Heisenberg spin— 1/2 chain in a zero 
magnetic field, 5*17(2) invariance fixes the exponents from purely phenomenological considerations. 



One of the central problems in condensed matter the- 
ory is the development of an effective phenomenologi- 
cal description of complicated many-body systems, the 
microscopic details of which are often not known. The 
description of low energy properties of interacting elec- 
trons in normal metals, for example, is provided by the 
theory of Fermi liquid [l| , while in one-dimensional sys- 
tems Luttinger liquid (LL) theory 0, |3| plays a sim- 
ilar role. These phenomenological theories do not rely 
on specific microscopic details, but predict certain low 
energy properties of many-body systems in terms of few 
measurable parameters. For example, LL theory assumes 
a linear spectrum of low energy excitations, and relates 
long-range behavior of correlation functions to the di- 
mensionless LL parameter K. However, the linear spec- 
trum approximation is not sufficient for finding dynamic 
response functions (DRFs) even in low energy limit [^. 
In this Letter we show, that for a wide class of ID systems 
one can phenomenologically predict certain properties of 
DRFs in terms of other measurable quantities even be- 
yond the low energy limit. 

The nonlinearity of the excitation spectrum affects 
transport phenomena, such as Coulomb drag Q and 
momentum-resolved tunneling of electrons Q between 
quantum wires. In addition, neutron scattering on 
spin chains ARPES on quasi- ID materials 0, and 
photoemission spectroscopy [10| of ID ultracold atomic 
gases directly measure DRFs, and are not limited to 
low energies. Evaluation of DRFs of ID quantum sys- 
tems with generic excitation spectrum is also a test 
bed for rapidly developing methods of numerical simu- 
lations of many-body dynamics [ll|, ■ Recently some 
progress was achieved in the analytical treatment of cor- 
relation functions beyond linear spectrum approxima- 
tion 0, [13, [3, [13, 0, [B, [13, Ei, [il H m ^ Hi]. 

The majority of analytical work however relied on solu- 
tions of microscopic models, using perturbation theory 
methods [13, [ll| or integr ability of models with specially 
tuned parameters [H [ll, [H, [H [ll, [11, [13, [HI . In con- 
trast, the phenomenology developed in this Letter does 
not require any special property of the underlying micro- 
scopic interaction, while it provides relations between dif- 



ferent experimentally observable quantities, such as the 
energy spectrum and the exponents of DRF singularities; 
see Eqs. (0) and (dH). 

For spinless fermionic Galilean-invariant systems, the 
DRFs have a sharp edge of support £{k) in the thermody- 
namic limit at T = 0; see Fig. [TJ Hamiltonian describing 
singularities of DRFs, such as dynamic structure factor 
(DSF) S{p, w) and spectral function A{p, uj) [defined be- 
low by Eqs. ([ij-lll])], is the effective Hamiltonian of a 
mobile impurity moving in a LL [l^, [B, [i3j [ISj [ll] ; see 
Eqs. ©-(O below. Singularities of DRFs at the edges 
of support are the main subject of this Letter. We show 
that their exponents for Galilean-invariant systems with 
interactions decaying faster than (xl/x can be expressed 
as functions of £{k) and LL parameters. Phenomenolog- 
ical considerations allow us also to resolve the discrep- 
ancy 12, T3| regarding antiferromagnetic spin— 1/2 XXZ 
model in zero magnetic field in favor of Ref. [l^ . 

We are interested in the zero temperature DSF 

S{k,Lo) = yda;die*("*^'=^' (p(a;,i)p(0,0)), (1) 



and spectral function A{k^uj) — — ilmG'(fcjw) signw, 
where Green's function G{k,uj) is defined by [2g| 



G{k,Lu) = -i 



da;(ite'("*-'^^)(r(*(a:,t)^'t(0,0))).(2) 



Here ^(a;,i) and p{x,t) are annihilation and density op- 
erators, respectively, and T denotes time ordering. En- 
ergy Lo is measured from chemical potential, so A{k,uj) 
for a; > (w < 0) describes the response of the system 
to an addition of an extra particle (hole). 

To be specific, we first discuss singularities of fermionic 
spectral function in the region |fc| < kf,uj < 0, where 
kf is Fermi momentum. Singularity can be described 
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FIG. 1: (Color online) (a) Spectral function A{k,ij) for spin- 
less fermions in momentum-energy plane at T = 0. Shaded 
areas indicate the regions where A{k,uj) ^ 0, and e(fc) is 
the edge of the support in the basic region |fc| < kf,uj < 0, 
where kf is Fermi momentum. Edges in other regions can 
be obtained from e{k) by combinations of shifts and inver- 
sions, (b) A sketch of constant k scan of the spectral function 
for |fc| < kf. Singularity of A(k,uj) near u ~ e{k) can be 
described using three-subband model of a mobile impurity 
moving in a Luttinger liquid, Eqs. Q-®, and the answer is 
given by Eqs. ©.(HIJ. 



13, 3 14 1 by the effective Hamiltonian 



Hf) = — I dx 
Zn 



Hd = J dx d){x) 



e{k) - ivd-Q^ 



dix), 



Hi, 



dx [Vrpr{x) + VLPLix)] pd{x) 
= J dx {Vn^'-^ - Vl^'-^) dix)dHx). 



(3) 
(4) 

(5) 



Here v is the sound velocity, and fields 6 and describe 
low energy excitations and have a commutation relation 
Wx),\/9{x')] = iwdlx — x') (we use the notation of Ref. 
[3|). Operator d{x) creates a mobile hole of momentum 
k and velocity Vd = de{k)/dk, and operator pdix) = 
d{x)d''{x) is the hole density. In terms of hole operator 
d{x,t), The singular part of the spectral function near 
e(fc) is given by 

A{k,uj)^ J dxdte'^*{d^{x,t)d{0,0))H„+H,+H,„,- (6) 

Canonical transformation — (f>\/K,d — diago- 



nalizes Hq, while the term Hint can be removed 25| by 
unitary transformation U'^ [Hq -\- Hd -\- Hint)U, where 



-ifdx{'-^ieix)-i(x)]- 



^-[4,{x)+0(x)]^d{x)d''{x) 



Momentum dependent phase shifts 6+{k), S-{k) are re- 
lated to the parameters of Hi„t as 



{Vl - Vr) /VK = -S^{k){vd + v) + 6+{k){va - v), (7) 
{Vl + Vr)^= -6-ik){vd + v)- S+{k){vd - v). (8) 

Calculating Wd{x)U together with Eq. ([6]), one obtains 
A{k,ULi) e{e{k) ~ u ^ 



e{k) — u) 



(9) 



To obtain phase shifts, one needs to fix Vr and Vl in 
Eq. (O. We relate Vr and Vl to e{k) by calculating 
in two ways the shift of the position of the edge under 
uniform density and current variations. 

Uniform density variation 6p results in a finite expec- 
tation value {Vlp) ~ —nSp. Evaluating the shift of e{k) 
in two ways, we obtain 

_ Vr + Vl _ de{k) dp _ de{k) wr 

2 ~ dp dp~ dp K' ^ ^ 

The left-hand side follows from the effective Hamiltonian 
given by Eqs. ©-(HI), while the right-hand side follows 
from the evaluation of the edge position from its thermo- 
dynamic definition [taking into account that energy e{k) 
is measured with respect to chemical potential p] . 

Uniform current through the system results in a finite 
value of (V0), which for Galilean-invariant systems corre- 
sponds to a motion with constant velocity u — {S/9) /m, 
where m is the bare mass of the constituent particles. 

27[, one can 



Then following the argument of Refs. [22 
use Galilean invariance to evaluate the change of e{k). 
Comparing it with the change evaluated using Eqs. ([3])- 
© leads to 



Vl - Vr _ _fc _ de{k) 
2tt m dk 



(11) 



Combining now Eqs. (O and ^ with Eqs. ITUl) and 
(fTTj) . we obtain the central result of this Letter 



5±{k) 
27r 



1 ( k_ 



dEjk) 
dk 



dejk) 
dp 



K 



(12) 



On the basis of Galilean invariance, it establishes a 
model-independent phenomcnological relation between 
the edge position e(fc) and other measurable quantities, 
such as the exponent of spectral function, Eq. ([9]). Even 
for usual LL theory Galilean-invariant systems are spe- 
cial. For them, LL parameter K can be expressed as 
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a renormalization of sound velocity compared to Fermi 
velocity Vf in the absence of interactions, K = Vf/v = 
■np/{mv). Our results are a generalization of the special 
properties of Galilean systems beyond low energy limit. 

While Eq. ^ doesn't hold on a lattice, Eq. (jTUl) still 
works, and will be used below for the XXZ model. One 
can formulate an analog of Eq. pT|) for LLs on lattices 
using the derivative of e{k) with respect to total flux 
through the system under periodic boundary conditions. 
Energies are easier to evaluate numerically than correla- 
tion functions, so our results can be used as a benchmark 
for numerical methods for evaluation of DRFs. 

Away from the basic region |fc| < kf^ix) < 0, positions 
of edges can be obtained from e(k) in the basic region 
by a combination of inversions and shifts. States which 
define the positions of the edges are given by a hole and 
excitations near Fermi points. Exponents of divergences 
can also be obtained using the three-subband model given 
by Eqs. ([3])-(l5]), and here we only summarize the results. 

For the spectral function momentum k in the region 
{2n — < k < (2n + l)fc/, hole momentum equals 

kn = k — 2nkf. Near the edges for uj > 0{llj < 0), the 
spectral function is defined by 



A{k,uj) oc 9{e{kn) ±uj) 



1 



LU ± £{kn) 



MTi.±(fc) 



27r 



1 ^1±1 , J +(fc„)-(5_(fc„) 
27r 



K 



(13) 



(14) 



DSF S{k,uj) is non-vanishing only for u; > 0, and for 
2nkf < k < 2{n + l)kf hole momentum is given by fc* = 
(2n + l)kf — k. The exponent of DSF /i„(fc) is defined by 



Sik,uj)oie{uj + e{k*)) 



{2n + l)VK- 



S+ik*J + s_{k;j 



K 



27r 

5+{kl)-5.{kl) 
2tt 



(15) 



.(16) 



For bosons, spectral function has divergences at 
=Fe(fc^) for Lo > 0{uj < 0), and exponents equal 



<±(fc) = l- o 



(2?i + 1)VK 

1 /1± 1 
'2 



5+{k:) + s^{k*„) 



K 



2tt 

S+{k*J-S.ik* 
2tt 



(17) 



Let us now discuss several cases where one can explic- 
itly check our phenomenological predictions. The shift of 
the position of e{k) can be evaluated using perturbation 
theory in interaction strength for any momenta, and pre- 
dictions following from our theory coincide with results of 



Refs. [l3, 14|. By using approximation e{k) « v{k — kf) 
for any interaction strength in the vicinity of the right 
Fermi point, from Eqs. jS]) and (fTOjl one can recover the 
universal phase shift [H] 



S-{kf-0) 
27r 



1 



2VK 



V' (18) 



which holds irrespective of Galilean invariance. One can 
also obtain (5+(fc/ — 0) from Eq. (fT^ . For that, one has to 
use the expansion e(fc) « i;(fc— fc/)-f (fc — fc/)f/(2TO*), and 
the expression for 1/m* obtained in Ref. [15|, which is 
valid for Hamiltonians with interactions decaying faster 
than oc 1/x^. For Galilean-invariant systems, it simpli- 
fies to 1/m* = [VK/{2'K)]dv/dn+l/{2mVK), and after 
some simple algebra with Eq. p2p one reproduces univer- 
sal phase shift (5+(A;/-0)/(27r) = 1-1/{2VK)-\^/2. 
One can also explicitly check, that exponents for Lieb- 
Liniger fiil and Calogero-Sutherland [2§| models evalu- 
ated from their excitation spectra reproduce the results 
of Refs. [3| and 14, 23|, respectively. 



The crucial step in the calculation of the exponents 
is the identification of the spectral function A(k,Lu) de- 
fined in terms of constituent particles, Eq. with the 
correlation function of operator d in Eq. Compari- 
son with solvable cases above shows that such identifica- 
tion indeed holds in the vicinity of Fermi points for any 
interactions, as well as for any momentum for weak in- 
teractions. While we cannot prove that it holds for any 
strongly interacting Galilean-invariant system, we expect 
it to be valid as long as the position of the edge satisfies 



d£{k) 



dk 



< V for |fc| < kf, 



(19) 



and interactions decay faster than cx l/x. Equation (jl9p 
guarantees that phases in Eq. (fT2|) are continuous func- 
tions of momentum, and the state which corresponds to 
the edge of the basic region of the spectral function sup- 
port does not contain particle-hole excitations near left 
or right Fermi points. 

Let us now discuss how considerations of this Letter 



can resolve a discrepancy between results of Refs. [12 1 
and 17| for correlations of the XXZ model in zero mag- 
netic field. In our notations, these references predict 
SPWA^f^y ^2Tr) = -5™(fc)/(27r) = {1/VK - ^)/2 
and 5^^(fc)/(27r) ^ -S^P{k)/{27r) = (1 - K)/2, respec- 
tively. Identification of these phase shifts was based on 
the analysis of finite size corrections to energies obtained 
from the exact solution. Their interpretation for the XXZ 
model in zero magnetic field is ambiguous, since the half- 
filled lattice is a special point for the Bethe Ansatz so- 
lution [s^l- On the other hand, our approach constrains 
phase shift via £{k), which is well defined in the thermo- 
dynamic limit, and resolves the discrepancy. 

First, we note that S^^^{k) satisfies universal relation 
given by Eq. (fT5|) . while 6^^{k) does not. Second, Eqs. 
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([7]), ([8]) and ([10]) hold on the lattice for any k, and one can 
easily evaluate excitation spectrum of the XXZ model nu- 
merically from the exact solution (30| . This way, we have 
verified that S^^^{k) satisfy them, while S^^{k) do not. 
Third, one can use SU (2) invariance to independently de- 
rive results of Ref. [l2| at the XXX point. The argument 
is very similar to the reasoningwhich fixes LL parameter 
K = 1/2 at XXX point 0,11] by requiring that long 
distance asymptotes of {S%x)S%0)) and {S-{x)S+{0)) 
coincide. But SU{2) symmetry also establishes a rela- 
tion between spin DRFs S^^{k,uj) and S ^{k,uj) in en- 
tire momentum-energy plane, including the edges of sup- 
ports. There S^^{k,uj) behaves as the sum of different 
power laws (up to logarithmic corrections) 



1 



to — e{k) 



A"! 



(20) 



Different power laws appear because of the umklapp pro- 
cesses that are allowed on a half-filled lattice. SU{2) 
symmetry implies that the same set of exponents should 
apply for S ^(k,Lu) as well. These exponents can be 
evaluated in terms of 6±{k) and K for any I, and the 
coincidence of two sets of exponents unambiguously fixes 

^=T^, (21) 

as in Ref. ^l?] for if = 1/2. Full sequence of exponents is 

;xf^ ^ 3/4 - (4; + 1) 2/4. (22) 

Note, that we did not use integrability in the argument 
for the XXX model, so Eq. ^ should apply for SU{2) 
invariant models with longer range interactions as well, 
if the spin chain remains a gapless LL. 

To summarize, we have considered zero temperature 
dynamic response functions of ID systems near edges 
of support in the momentum-energy plane. Continuous 
symmetries can be used to fix the exponents of power law 
divergences of dynamic response functions near the edges. 
For spinless Galilean-invariant systems of fermions or 
bosons, we have obtained phenomenological expressions, 
Eqs. (in]),([Tl]),([ini), and which establish model- 

independent relations of the exponents of dynamic re- 
sponse functions to the position of the edge of support 
e(fc). For a spin— 1/2 anitferromagnetic Heisenberg chain 
in zero magnetic field, SU{2) symmetry dictates expo- 
nents given by Eq. (|^^ for all momenta regardless of the 
interaction range. 
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